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NONLOCAL  EFFECTS  IN  BRITTLE  CRACK  PROPAGATION* 


N.  Ari  and  A.C.  Eringen 
Princeton  University 
Princeton,  NO  08544 


ABSTRACT 

Stress  distribution  near  the  tip  of  a  constant  velocity  crack  is 
determined  by  means  of  nonlocal  theory  of  elasticity.  The  maximum  stress 
is  finite  and  it  allows  one  to  utilize  the  maximum  stress  hypothesis  to 
determine  the  terminal  velocity  of  running  cracks. 


1.  INTRODUCTION 

One  of  the  fundamental  mechanisms  of  dynamic  rupture  phenomena  is 
that  of  crack  propagation.  The  angular  distribution  and  the  maxima  of  the 
crack  tip  stresses  vary  strongly  with  the  crack  velocity  and  the  dynamic 
effects  become  important.  One  has  to  solve  an  elastodynamic  problem  to 
obtain  these  dynamic  stresses  so  that  at  least  a  qualitative  understanding 
of  the  crack  motion  can  be  understood. 

The  present  work  investigates  stress  distribution  near  the  tip  of 
a  moving  crack  in  a  brittle  elastic  solid  by  means  of  the  recently  developed 
theory  of  nonlocal  elasticity  [1]. 
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The  dynamic  crack  problems  are  rather  difficult  to  solve  without 
some  simplifying  assumptions.  Most  common  among  them  are  the  constant 
crack  tip  velocity  and/or  the  loading  conditions  which  allow  self-similar 
crack  motions.  As  such,  these  problems  represent  a  highly  contrived  picture 
of  the  actual  rupture  phenomenon.  Nevertheless,  these  models  prove  to  be 
useful  in  discussing  the  velocity  dependence  of  the  crack  tip  stresses. 

One  of  the  earliest  and  simplest  models  of  a  moving  crack  is  the 
Yoffe's  model  [2].  A  straight  line  crack  of  fixed  length  21  moves  with 
a  constant  velocity  in  an  infinite  plate  subject  to  uniform  tension  perpen¬ 
dicular  to  the  crack  (Fig.  1).  Below,  we  present  a  solution  of  Yoffe's 
problem  within  the  context  of  nonlocal  elasticity. 

The  classical  elasticity  solution  of  Yoffe's  problem  as  well  as 
other  classical  treatments  of  uniformly  moving  cracks  (3-5)  have  provided 
certain  insights  into  dynamical  aspects  of  brittle  fracture  and  inspired 
subsequent  research.  However,  these  studies  have  remained  inadequate  in 
explaining  a  number  of  important  features  of  dynamic  crack  propagation. 

To  this  end,  the  following  points  are  worth  noting: 

(i)  The  classical  elasticity  solution  yield  the  familiar  square  root 
stress  singularities  at  crack  tips.  In  order  to  circumvent  the 
difficulties  related  to  infinite  stresses,  the  dynamic  fracture 
criteria  are  based  on  energy  considerations  (energy  release  rate). 
However,  the  dynamic  failure  is  localized  at  the  crack  tip  and 


it  depends  on  the  critical  stress  levels.  A  maximum  stress  criterion 


is  therefore  more  appropriate. 
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(ii)  The  dynamic  stress  intensity  factor  in  Yoffe's  solution 

(1.1)  K  =  Lim  [(27rr)^  t  (r,e=0)] 

r+0  yy 

turned  out  to  be  independent  of  crack  tip  velocity  and  equal  to 
its  static  value.  This  result  is  contrary  to  experimental  obser¬ 
vations  that  K  should  vary  with  the  crack  speed. 

(iii)  The  analytic  expressions  of  the  classical  solutions  imply  that 
the  Rayleigh  wave  velocity  cR  is  a  natural  upper  bound  for  the 
crack  tip  velocity.  However,  for  all  materials  observed,  the 
terminal  velocity  V^.  remained  considerably  less  than  cR  [6], 

Since  Mott  [7]  first  addressed  the  problem  of  a  modified  Griffith  • 
criterion  for  a  moving  crack,  the  determination  of  the  limiting 
velocity  Vt  has  been  an  open  problem.  Operating  on  the  premise 
of  a  constant  surface  energy  Yq  ,  Roberts  and  Wells  arrived  at  a 
formula  expressing  the  terminal  velocity  [8],  However,  experimental 
evidence  indicates  that  the  surface  energy  varies  with  the  crack 
velocity,  i.e.  y  =  y(V).  Therefore,  Freund  [9]  considers  any 
agreement  of  their  theoretical  terminal  velocity  with  the  experimental 
one  rather  accidental.  An  additional  difficulty  in  dealing  with  sur¬ 
face  energy  concept  is  its  experimental  determination  and  accuracy. 

(iv)  A  more  involved  fracture  phenomenon  is  that  of  crack  bifurcation  — 
a  rapidly  moving  crack  suddenly  branches  into  two  new  cracks.  The 
precise  physical  source  of  this  mechanism  is  not  known.  Neither  is 
it  settled  as  to  which  criterion  is  better  suited  to  determine  the 


prebranching  velocity  and  the  bifurcation  angles  once  the 
solution  for  the  dynamic  problem  is  given  [10]. 

The  motivation  for  the  present  work  is  derived  from  the  foregoing 
considerations.  In  order  to  develop  controlled  fracture  processes  and  crack 
arrest  mechanisms,  the  need  for  addressing  the  above  questions  is  clear.  En¬ 
couraged  by  the  results  of  the  static  nonlocal  crack  problems  [11-13],  we  will 
pursue  the  solution  of  the  nonlocal  Yoffe  problem.  The  nonlocal  theories 
incorporate  the  long  range  interactions  and  the  microstructure  dependence 
intrinsically.  Hence,  they  provide  a  more  realistic  framework  for  the 
treatment  of  crack  problems. 

Finally,  the  formulation  of  a  mixed  boundary  value  problem  in  nonlocal 
elastodynamics  is  in  itself  important.  Difficulties  related  to  mixed  boundary 
conditions  have  been  considered  in  our  earlier  work  [12,14],  In  the  present 
work,  in  addition,  we  discuss  a  modification  of  the  field  equations  due  to 
the  dynamic  effects.  Throughout  the  work,  emphasis  will  be  on  how  the 
nonlocal ity  modifies  the  dynamic  crack  tip  stresses. 

Section  2  formulates  the  nonlocal  boundary  value  problem.  In  Section 
3,  the  ensuing  dual  integral  equations  are  solved.  The  velocity  dependence 
of  the  crack  tip  stresses  is  given  in  Section  4.  Section  5  compares  the 
nonlocal  results  with  the  classical  ones  and  discusses  their  relevance  to 
dynamic  fracture  criteria.  The  quantitative  results  are  gratifying. 

2.  NONLOCAL  FORMULATION 

In  nonlocal  elasticity  for  isotropic  solids,  the  stress  constitutive 
equations  are  given  by  [1]. 
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(2.1) 


[X * ( | x ' -x | )  err(x-)£kS,  +  2jj * ( | x ' -x j  )ek£ (x* )]da(x' ) 

.  J 


where  the  integral  is  over  the  two-dimensional  plane  region  R(x',y')  and 
e.  is  the  strain  tensor  which  is  related  to  the  displacement  vector  uk  by 


(2.2) 


6ki"  I  K,£  +  U£.k} 


Here,  and  throughout  an  index  following  a  comma  represent  gradient, 


uk,«.  ■  3V3xt 


The  nonlocal  kernels  X * ( | x * -x j )  and  u'(|x'-x[)  represent  the 

*** 

influence  of  the  neighboring  strains  at  x‘  on  the  nonlocal  stress  at  a 
reference  point  x.  They  are  usually  determined  by  requiring  that  the 
nonlocal  field  equations  yield  identical  dispersion  relations  to  those 
derived  from  atomic  lattice  dynamics.  In  this  way,  in  [12],  we  introduced 
a  two-dimensional  kernel 


(2.3)  X'/X  =  y'/u  =  °t(  | x '  —x | )  =  ^  62  K0^[(x,-x)2  +  (y'-y)2]^> 


where  Kg  is  the  zeroth-order  modified  Bessel  function  and  6  is  the 


nonlocality  parameter.  In  the  local  limit 


(2.3)  reverts  into 


a  Dirac-delta  functional  and  (2.1)  yields  the  classical  Hooke's 
law.  In  addition,  (2.3)  has  the  convenient  property  that 
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(2.4)  (1  -  e2V2)  cx(x,y)  =  <5(x,y) 

where  e  =  l/g  and  <5  is  the  Dirac-delta  functional. 

For  nonlocal  elasticity  (when  the  effects  of  nonlocal  residuals  and 
the  body  forces  are  excluded)  the  integral  balance  equations  of  linear 
momentum  yield  the  equations  of  motion 

(2*5^  tk£,k  =  pu£ 

where  a  superposed  dot  represents  the  material  time  derivative.  Due  to  the 
integral  form  of  constitutive  equations,  (2.5)  is  of  integro-differential 
character.  For  such  systems,  it  is  very  difficult  to  establish  existence  and 
uniqueness  theorems  and  special  care  is  needed  in  the  formulation  of  mixed 
boundary  value  problems.  In  addition,  a  closer  look  at  the  field  equations 
signals  certain  difficulties  for  cases  where  the  displacement  fields  are  dis¬ 
continuous  such  as  in  bodies  with  cracks.  For  example,  in  (2.5),  when 
term  is  discontinuous,  the  t^  k  term  is  very  likely  not  to  be  discontinuous, 
since  the  discontinuities  in  the  displacement  fields  will  be  smoothed  out  by 
the  integral  operator  in  (2.1). 

With  the  appearance  of  the  crack,  some  interatomic  bonds  are  eliminated 
and  inhomogeneities  arise  within  a  narrow  boundary  layer  near  the  crack  sur¬ 
faces.  The  nonlocal  body  forces  become  important  and  the  use  of  nonhomogeneous 
kernels  are  required.  Since,  however,  the  inhomogeneities  are  confined  to  a 
narrow  domain,  we  can  approximate  the  problem  with  a  homogeneous  kernel,  pro¬ 
vided  the  compatibility  of  the  crack  discontinuities  with  the  nonlocal  field 
equations  are  ensured. 
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To  this  end,  a  modification  of  the  field  equations  can  be  obtained  by 
utilizing  the  relation  (2.4),  i,e, 

(2.6)  U-eV)  tkt>k  =  ok),jk  =  (1-eVl  p'J, 

where 

°kn  =  Aerr  6k«,  +  2u  ek£ 

is  the  local  (classical)  stress  tensor.  (2.6)  together  with  the  appropriate 
side  conditions  constitute  a  sufficient  solution  for  the  integral  balance 
equations.  The  modified  equations  of  motion  are  singularly  perturbed  dif¬ 
ferential  equations  and  they  are  simpler  than  the  integro-differential  system 
given  by  (2.5).  (2.6)  can  be  derived  alternatively  from  integral  balance 

equations  by  imposing  uniform  continuity  requirements  on  nonlocal  stresses 
and  displacements  [15]. 

In  the  two-dimensional  formulation,  we  assume  that  the  displacement 
fields  can  be  derived  from  wave  potentials  <f>(X,Y,t)  and  ip(X,Y,t),  i.e. 

(2.7)  u  =  dp/dX  +  3ij;/9Y 

v  =  dp/dY  -  dp/dX 

Substituting  (2.7)  into  (2.6),  we  obtain  two  scalar  wave  equations 

(2.8)  (1  _  e2v2)  <p  =  c2V2  p 

(1  -  e2V2)  $  -  CjV2  ip 

where 

c]  =  C( A  +  2y)/pT*  ,  c2  =  (p/ qY* 


mammssm 


which  must  be  solved  to  determine  the  displacement  field.  In  the  case 
of  the  constant  velocity  crack  problem,  it  is  more  convenient  to  intro¬ 
duce  a  moving  coordinate  system  (Fig.  1), 
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(2.9)  x  =  X  -  Vt  y  «  Y 

where  V  Is  the  crack  tip  velocity.  We  assume  that  in  the  new  coordi¬ 
nate  system,  the  wave  potentials  become  independent  of  time.  Then  the 
field  equations  (2.8)  reduce  to 


(2.10) 


3x  3x^  3y^  3x  3y^  3x^ 


<t>(x,y) 


c2V2  ^7  (^y  +  ^j)  *  c2(^  +  -  V2 

3x*  3x‘  3y^  c  3x^  3y^  3x^ 


*(x,y) 


0 


0 


By  means  of  the  Fourier  transform  in  the  x-direction,  we  obtain 
the  general  solution  of  (2.10)  satisfying  the  conditions  0  as 
(x2+yV/2+~ 


<J>(k,y)  *  A(k)  e 


"Y^(k)y 


(2.11) 


-Y2(*Oy 


<P(k,y)  *  B(k)  e 
where  k  is  the  transform  variable  and 


y  >  0 


(2.12) 


r-(k)  *  [U- -V2-e2k2V2)/(c2-c2k2V2)]k2  ; 


i*l  ,2 
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The  inversion  contour  of  (2.11)  requires  proper  branch  cuts  so  as  to  ensure 
that 

(2.13)  Re  y1 (k)  >  0 

The  boundary  conditions  for  moving  crack  in  the  new  coordinate 
system  is  equivalent  to  the  static  case.  Therefore,  the  self-consistent 
nonlocal  boundary  conditions  are  given  by, [12]  , 


(2.14) 


oyy(x,0) 
v(x ,0) 

axy(x.O) 

u,v  -*-0 


-  t0  |x|  <  £ 

0  |  x  |  >  £ 

0  V  x 

as  (x2+y2)^  -*■  ® 


(2.14)  and  (2.11)  lead  to  a  set  of  dual  integral  equations  whose  solution 
will  determine  A(k)  and  B(k). 


3,  SOLUTION  OF  THE  DUAL  INTEGRAL  EQUATIONS 


Substituting  (2.11)  and  (2.7)  into  (2.14)*  we  obtain 


(3.1)  oyy(x,0) 


oo 


{[(A  +2u)  Y]  -  Ak2]  A(k)  -  2U  ik  y2  B(k)}  e"kxdk 

-00 


=  -  t0  ;  i  X I  <  £ 


(3.2)  v(x,0)  =  (2tt) 


[-Y1  A(k)  +  ik  B(k) ]  e-lkx  dk  =  0  ; 

lx]  >  £ 


(3.3)  axy(x,0)  =  (2tt)" 


J 

-oo 


u(2ik  y1A(k)  +  (y|  +k2)  B(k)  e'lkx  dk  =  0 

;  V  x 


An  exact  solution  to  the  system  is  yet  to  be  found.  The  problem  is  further 
complicated  by  the  fact  that  the  Y^(k)  of  (2J2)  assumes  complex  values 
along  the  real  axis  between  their  respective  branch  points. 

In  the  static  crack  problems,  it  has  been  shown  that  the  classical 
solution  constitutes  a  reasonable  approximation  to  the  nonlocal  one  [11,12]. 
The  only  difference  between  the  classical  and  nonlocal  versions  of  the  system 
(3.1) -(3.3)  is  the  dependence  of  y..(k;  e)  to  e  ,  where  y^k,  c=0)  * 

(3.4)  s2  =  1  -  (V/Ci)2  * 

In  the  limit  as  e  -*■  0  ,  we  recover  the  classical  solution.  The  behaviour 
of  integrands  in  (3.1) -(3.3)  for  large  k  is  similar  to  the  classical  case, 
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therefore,  it  is  expected  that  the  limiting  process,  as  e  -*■  0  ,  shall 
be  uniform.  Hence  as  an  approximation,  one  might  use  the  classical  solution 
provided  that  the  complex  values  of  y^(k)  are  taken  into  account.  To  this 
end,  we  choose  a  modified  form  of  the  classical  solution  such  that 

( k)  AN(k)  and  y^(k)  B^(k)  remain  always  real  valued  and  numerically 
equal  to  the  absolute  values  of  y^  and  y^Bc  respectively.  These 
conditions  will  be  satisfied  if  the  nonlocal  solutions  for  AN  and  BN  are 
expressed  as 


(3.5) 

AN(k)  =  P^k)  +  1  P2(k) 

(3.6) 

BN(k)  =  Q-j  ( k )  +  i  Q2(k)  , 

where 

(3.7) 

P1(k) 

=  Ac(k)  Q-| (k)  =  Bc(k)  ; 

(3.8) 

P2(k) 

=  -  Ac(k)  Q2(k)  =  -  Bc(k); 

and  R  correspond  to  those  regions  of  the  real  k-line  on  which  y-j(k)  or 
Y2(k)  are  not  real  valued. 

The  classical  solutions  are  given  by  [16] 


(3.9) 

Ac(k)  = 

(1  +  s2)  D^kjO/k2 

(3.10) 

Bc(k)  = 

-2s,i  DpO^kJO/k2 

where 
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Dq  =  (f)^  t0  ^(V/C2)2  [4Sls2  -  (l-s2)2]’1 

It  can  be  shown  by  substitution  that  the  nonlocal  boundary  conditions  are 
uniformly  approximated  by  (3.5) -(3.6)  to  the  order  of  0  (e^)  [15]. 

In  the  present  work,  we  are  primarily  interested  in  the  normal  stress 
along  the  crack  line,  thus  in  the  next  section,  we  utilize  this  modified 
classical  solution  to  express  the  normal  stress  along  the  crack  line  near 
the  tip. 


4.  NONLOCAL  STRESSES 

The  constitutive  equations  for  the  stress  field  are  given  by  (2.1). 
Using  the  standard  integration  formulas  and  substituting  ( 3 . 5) -( 3. 1 0)  into 
(2.1),  we  obtain  for  the  normal  stress  along  the  crack  line 


(4.1) 
where 

(4.2) 


(4.3) 


(4.4) 

and 


yy 


k-j [4  s  1  tg  -  2  (l+s‘)  tA] 


{tA’V  = 


0-j ( k)  cosk£  (aA(ck),  aB(ek)}  dk 


(C, -n2V2  -  e2k2V2)  k  (C2-  V2-  e2k2V2 

aA(ek)  =  -1 - ^ - [1  -  TT^TJTk  d  "72 - TTJrO  3 


'1 


(l+e2k2)^  (cf  -  e2k2V2) 


aB(ek) 


C2  -  V2-e2k2V2 
{(| - 5 - T~T~5 - 1)  "  E 


C2  -  e2k2v2 


k  C2  -  V2  -  e2k2V2  C2  -  e2k2V2 

TO  3 - 


(l+eV)'1  C~  -  e  k  V 


2\2t  j.  n  %/ 


kl  =  E4  -  (1  +  s2)‘]  tQ  2tt 
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(4.5)  n2  a  (A+2y)/2y 

C  =  x/£  ;  e  =  1/32. 

Near  the  crack  tip,  stresses  are  given  in  terras  of  generalized  hypergeometric 
functions  which  are  difficult  to  visualize  in  graphic  terms  (see  Appendix  1). 
In  the  present  work,  we  are  most  interested  in  the  velocity  dependence  of 
the  crack  tip  stresses.  Therefore  we  present  the  numerical  values  of 
tyy(l,0;  V)  in  Table  1. 

An  alternative  method  of  solution  to  the  system  (3.1)  -  3.3)  is  the 
reduction  of  the  dual  integral  equations  into  an  equivalent  Fredholm  in¬ 
tegral  equation  of  the  second  kind.  The  Fredholm  integral  was  then  solved 
numerically.  The  crack  tip  stresses  so  obtained  are  comparable  to  the 
analytic  results  as  shown  in  Table  2. 
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5.  DISCUSSION 

The  following  results  are  worth  noting: 

(i)  The  nonlocal  stresses  at  the  crack  tip  are  finite.  These  results 
will  enable  us  to  extend  Eringen's  fracture  criterion  [17]  to 
dynamic  rupture  phenomena.  In  the  limit  as  the  crack  velocity 

V-"  0  ,  we  recover  the  static  nonlocal  crack  results  (see  Table  1). 

(ii)  The  local  analysis  of  the  steady-state  crack  problems  yielded  stress 
intensity  factors  which  are  independent  of  the  crack  velocity  and 
equal  to  their  static  values  (Yoffe  [2]  Craggs  [3]).  These  results 
are  clearly  undesirable  in  that  they  contradict  the  experimental 
results  which  show  that  the  crack  tip  velocity  depends  on  applied  loads. 

The  nonlocal  analysis  on  the  other  hand,  provided  a  crack  tip  stress 
behaviour  which  is  compatible  with  experiments.  The  substantial 
difference  between  the  local  and  nonlocal  results  can  be  observed 
in  Figure  2. 

(iii)  The  nonlocal  results  can  be  employed  further  to  determine  the 
terminal  velocity  of  propagating  cracks.  The  uniformly  moving 
crack  is  in  a  state  of  dynamic  equilibrium.  As  the  crack  pro¬ 
pagates  first,  the  interatomic  bonds  at  the  tip  break.  Then 
as  the  stress  around  the  crack  tip  relaxes,  the  next  bond  is 
overloaded  by  the  sweeping  stress  wave.  There  will  be  a  sudden 
increase  in  the  stress  level  above  the  ultimate  stress  one 
atomic  distance  away  from  the  tip. 
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At  the  same  time,  from  Fig.  2,  we  observe  that  the  crack  tip 
stress  decreases  with  increasing  velocity.  In  the  dynamic  equi¬ 
librium,  two  changes  in  the  stress  level  compensate  for  each 
other  and  the  crack  reaches  a  uniform  propagation  velocity. 


The  stress  relaxation  at  the  crack  tip  is  a  complex  phenomenon. 
Therefore,  it  is  difficult  to  calculate  the  stress  overshoot 
precisely.  However,  using  a  harmonic  oscillator  model,  we  may 
assume  that  it  is  roughly  equal  to  the  difference  between  the 
crack  tip  stress  and  the  stress  level  one  atomic  distance  away 
from  the  tip.  Consequently,  we  propose  that  the  terminal 
velocity  be  determined  by 


(5.1) 


t¥¥(l;Vt) 

y> ;(»  '  t"r. ;  'o 


1S  given  in  £12,  Table  1]  for  the  static  crack 
problem.  tyy(v)/tyy(°)  is  plotted  on  Fig.  2.  On  Fig.  3,  both 
curves  are  drawn  together  —  their  crossing  point  yields  the  value 
of  the  terminal  velocity  as 


(5.2) 


Vt  -  pc2  pM0.5  “  0.6) 


Note  that,  in  the  case  one  is  able  to  solve  the  stress  relaxation 
problem  at  the  crack  tip,  ^yyd^J/tyyO)]  ™  fl‘9*  3,  may  not 
be  a  straight  line  any  more  (i.e.  it  will  be  velocity  dependent). 
In  that  case,  the  terminal  velocity  estimates  will  be  slightly 
modified. 


(iv)  By  formal  extension  of  Eringen's  fracture  criterion  to  dynamic 
phenomena,  we  can  predict  qualitatively  the  variation  of  the 
surface  energy  with  the  propagation  velocity.  The  Griffith 
definition  of  surface  energy  can  be  formally  extended  to  the 
velocity  dependent  case  by 

(5.3)  t\  a  =  -^jy  C(V)2  Y (V)  , 

where  y(v)  is  the  surface  energy,  p,v  are  the  standard 
material  moduli,  tj.  is  the  cohesive  stress,  a  is  the  inter¬ 
atomic  distance  and  C(V)  is  defined  by 

(5.4)  tyy(l,0)/t0  *  C ( V)  (2£/a)V2  . 

By  assuming  the  cohesive  strength,  a  material  property  and 
therefore  constant,  from  (5.3)  we  can  deduce  that 

<«■»  S  *  C§[2|32 

From  (5.3)  and  Table  1,  we  can  deduce  that  the  surface 
energy  increases  with  increasing  crack  velocity.  Although  there  is 
considerable  scatter  in  the  experimental  values  of  dynamic  surface 
energies,  the  nonlocal  results  are  in  agreement  with  their  general 
tendency  to  increase  with  crack  speed  [9], 
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(v)  If  one  wishes  to  establish  a  link  with  classical  fracture  theories, 
the  following  is  interesting  to  mention.  A  typical  fracture  tough¬ 
ness  versus  crack  velocity  relationship  is  shown  in  Figure  4  [18]. 

The  general  character  of  this  relationship  is  supported  by  physical 
considerations  (Erdogan  [19],  Freund  [20]),  as  well  as  by  experiments 
on  structural  mild  steel  [21,  22]  and  on  glassy  polymers  [23].  On 

Figure  4,  we  also  plotted  the  ratio  t  (V)/t  (0). 

^max  yymax 


Figure  4  indicates  that  in  the  lower  velocity  ranges,  the  fracture 
toughness  decreases  faster  than  the  maximums  tress  with  increasing  velocity. 
Thus  when 


(5.6) 


yy, 


(v)  >  t  (v) 

max  critical 


» 


the  crack  will  accelerate.  According  to  Eringen's  fracture  criterion,  a 
dynamic  equilibrium  will  be  established  when 


(5.6) 


W'*’  ■  ^critic, 


We  can  therefore  determine  the  terminal  velocity,  assuming  that  it  is 
defined  by  this  dynamic  equilibrium.  From  Figure  4,  we  observe  that  the 
terminal  velocity  will  be  on  the  order  of 


(5.7) 


Vt  n,  0.55  C2  . 
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The  important  feature  of  this  result  is  in  that  it  predicts  a  terminal 
velocity  well  below  the  Rayleigh  wave  speed, which  is  consistent  with 
experimental  observations  and  some  numerical  simulation  studies  [24]. 
Furthermore,  it  is  obtained  by  a  simple,  natural  fracture  criterion, 
without  extraneous  arguments.  The  present  approach,  however,  should  be 
considered  with  some  caution,  since  the  fracture  toughness  measurements 
include  some  global  effects  such  as  viscoelasticity  and  plasticity.  There¬ 
fore,  we  prefer  the  microstructural  considerations  discussed  in  (iii), 
to  predict  the  terminal  velocity. 

(vi)  The  explanation  of  crack  bifurcation  phenomena  is  much  more 
difficult.  Experimental  observations  suggest  that  [25] 

(1)  °f  Cb  =  constant 

where  af  is  the  fracture  stress  and  Cb  is  thelength  of 
the  crack  at  the  moment  of  branching. 

(2)  "Crack  branching  is  not  a  spontaneous  event"...  but  is 

the  eventual  outcome  of  cumulative  process  of  advance 
cracking. 

These  characteristics  of  crack  branching  suggest  that  it  is  not  likely 
to  explain  crack  bifurcation  solely  in  terms  of  a  critical  branching 
velocity.  One  has  to  consider  the  fracture  energies  of  the  involved 
materials  (i.e.  glasses  do  have  low  but  tool  steels  do  have  high  fracture 
energies).  In  addition,  the  accelerating  phase  of  the  crack  propagation 
has  to  be  taken  into  account. 

Vor  example,  Broberg's  local  analysis  [4]  predicts  Vt  =  V^. 
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Our  present  analysis  does  not  include  considerations  about  advance 
cracking  and  acceleration  effects.  Hence  we  do  not  attempt  to  explain 
crack  branching  on  the  basis  of  the  present  result.  Nevertheless,  from 
the  foregoing  conclusions.  It  is  clear  that  the  nonlocal  dynamic  results 
are  in  the  right  direction. 


TABLE  IV.  1 

VARIATION  OF  CRACK  TIP  STRESS 
WITH  THE  CRACK  TIP  VELOCITY 

(from  4.3,  X  ) 


t  (l,0;V)/t  (1,0;  V*0) 


I  I  t  t  t  r 


FIGURE  1 

PLATE  WITH  MOVING  LINE  CRACK 


Classical  Solution  [2] 


APPENDIX  1 


The  evaluation  of  the  stress  tip  integrals  can  be  facilitated  by 
the  use  of  the  asymptotic  expansion  of  Bessel  functions,  i.e. 


(A.l) 


J-j { k)  *  k_%  cos  (k  -  3tt/4)  +  0(k'3/2)  ; 


Substituting  (A.l)  into  (4.2),  we  obtain  the  following  results  as  £-*•  1 


(A. 2)  tk  *  C(cf  -  n2v2)/C2]  (ta  +  tc  +  te)  -  (tb  ty  tc) 


(A.3)  tB  =  t  ♦  t,  +  t, 


where 


t  7  n  *5/2  M3/2  2——  r(3/2)  r(3/4) 

(A. 4)  ta  -  2  s]  M1  -  0.5  s1  M1  5/4  r(g/8)  r(13/8) 


3F2(l/2,  3/2,  3/4;  13/8,  9/8;  qa) 


+  _  n  ,  5/2  J/2  n  ,  9/2  m3/2  tl_  r(3/2)  r(7/4 

(A. 5)  tb  -  0.4  s1  M]  -  0.5  s1  M1  9/4  r(l3/8)  r(l7 


-s|  M2/ 4 


"1  29/4  r(13/8)  r(17/8) 

3F2  (1/2,  3/2,  7/4;  17/8,  13/8;  q#) 


(2.6)  tc  =  2(l-sJ/2)  M1/2  -  0.5  M^3/2  sj1/2  B(l/2,  3/2) 


I  £' 

n=0 


(l/2)„  (3/2) 


-  Mj2n  ^(1/4,  n+3/2;  n+2;  qc)] 
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(A. 7)  t 


(A.8)  t, 


(A. 9)  t 


(A.  10)  t 


(A. 11)  t 


d  =  0.4  (1-s^2)  m]/2  -  0.5  M‘3/2  s3/2  B ( 1 / 2 »  3/2) 

0/2)  (3/2) 


l  c- 

n=0 


%  *  '  Hi2  si' 


'n 

T27 


n 


-  M‘2n  2F1(-3/4,n+3/2;n+2;qc)] 


0.5  ±S!i"  (1/2)"  M/4)» 


TT74V 


2F,(-l/2,  n-1/4;  n+1/4;  s2)]] 


0  5  H,/2  TO  1  t  r<!/2!  F(-5/4)  r  y  tl/2)n  (_5/4)n  (’H12) 
°,5  H,  [0.4  ♦  ■  [Jo - P37477- 


n! 


2^i (-1/2,  n-5/4;  n-3/4;  s2)] 


=  0.5  M^/2[s3/2  B( 3/2 ,  1/4)  2F1 (-1/2,  1/4;  7/4;  s2) 

-  s7/2  M2  B(2 ,3/4)  2F1(l/2,  3/4;  11/4;  -  s|/m|)] 


=  0.5[M27/2  s23/2  B ( 3/ 2,3/2)  ^(3/4, 3/2;  3;  -  (M^)'2 

-  [(4/7)  M"2  q7/4  F(7/4,9/4;ll/4;qi)  -  s7/2  ^ ( 7/4 ,9/4 ;1 1 /4 ;s|) 

-  (4/3)  s2  q3/4  F ( 3/4 , 5/4 ; 7/4 ;q. )  -  s3/2  2F1(3/4,5/4;9/4;s2)]] 

-•  C2/(C2  +  V2) 
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Air  Dniversity  Library 
Maxwell  Air  Forca  Base 
Montgomery,  Alabama  36112 

Other  Government  Activities 

Commandant 

Chiaf,  Tasting  and  Development  Diviaion 
8.S.  Coast  Guard 
1300  E  Street,  HW. 

Washington,  D.C.  20226 


Technical  Director 
Marine  Corps  Development 
and  Education  Command 
Quantico,  Virginia  22134 

Director  Defense  Research 
and  Engineering 
Technical  Library 
Room  3C128 
The  Pentagon 
Washington,  D.C.  20301 

Dr.  M.  Gaus 

Rational  Science  Foundation 
Environmental  Research  Diviaion 
Washington,  D.C.  20550 

Library  of  Congress 

Science  and  Technology  Division 

Washington,  D.C.  20540 

Director 

Defense  Mac leer  Agency 
Washington,  D.C.  20305 
Attn:  SPSS 
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Other  Government  Activities  (Con't) 

Mr.  Jerome  Penh 
Staff  Spaeialiat  for  Matariala 
and  Structures 
OUSDR&E ,  Tha  Pant agon 
Room  301089 

Washington,  D.C.  20301 

Chiaf,  Airframe  and  Equipment  Branch 
FS-120 

Offica  of  Flight  Standards 
Fadaral  Aviation  Agency 
Washington,  D.C.  20S33 

National  Academy  of  Sciaticas 
National  Sasaareh  Council 
Ship  Hull  Research  Committee 
2101  Constitution  Avenue 
Washington,  D.C.  20418 
Attn:  Mr.  A.  L  Lytle 

National  Science  Foundation 
Engineering  Mechanics  Section 
Division  of  Engineering 
Washington,  D.C.  20330 

Picatinny  Arsenal 

Plastics  Technical  Evaluation  Canter 
Attn:  Technical  Information  Section 
Dover,  New  Jersey  07801 

Maritime  Administration 
Office  of  Maritime  Technology 
14th  and  Constitution  Avenue,  NW. 
Washington,  D.C.  20230 

PAST  2  -  Contractors  and  Other  Technical 
Collaborators 

Pnivsreitiae 

Or.  J.  Tinsley  Odea 
University  of  Texas  at  Austin 
343  Engineering  Science  Building 
Austin,  Texas  78712 

Professor  Julius  Miklomitx 
California  Institute  of  Technology 
Division  of  Engineering 
and  Applied  Sciences 
Pasadena,  California  91109 


Universities  (Con't) 

Dr.  Harold  Liebovitz,  Dean 
School  of  Engineering  and 
Applied  Science 
George  Washington  University 
Washington,  D.C.  200S2 

Professor  Eli  Sternberg 
California  Institute  of  Technology 
Division  of  Engineering  and 
Applied  Sciences 
Pasadena,  California  91109 

Professor  Paul  M.  Naghdi 
University  of  California 
Department  of  Mechanical  Engineering 
Berkeley,  California  94720 

Professor  A.  J.  Durelli 
Oakland  University 
School  of  Engineering 
Rochester,  Missouri  48063 

Professor  F.  L.  DiMaggio 
Columbia  University 
Department  of  Civil  Engineering 
New  fork,  New  York  10027 

Professor  Norman  Jones 

The  University  of  Liverpool 

Department  of  Mechanical  Engineering 

P.  0.  Box  147 

Brown low  Hill 

Liverpool  L69  3BZ 

England 

Professor  E.  J.  Skudrzyk 
Pennsylvania  Stats  University 
Applied  Research  Laboratory 
Department  of  Physice 
State  College,  Pennsylvania  16801 

Professor  J.  Klosner 
Polytechnic  Institute  of  New  York 
Department  of  Mechanical  and 
Aerospace  Engineering 
333  Jay  Street 
Brooklyn,  New  York  11201 

Professor  R.  A.  Schapery 
Texas  A&M  University 
Department  of  Civil  Engineering 
College  Station,  Texas  77843 


Universities  (Con't.) 

Professor  Walter  0.  Pilkey 
University  of  Virginia 
Research  Laboratories  for  the 
.  Engineering  Sciences  and 
Applied  Sciences 
Charlottesville,  Virginia  22901 

Professor  t.  0*  Villmert 
Clarkson  College  of  Technology 
Department  of  Mechanical  Engineering 
Potsdam,  Hew  York  13676 

Dr.  Walter  E.  Haisler 
Texas  A&M  University 
Aerospace  Engineering  Department 
College  Station,  Texas  77843 

Dr.  Sustain  A.  Kamel 
University  of  Arizona 
Department  of  Aerospace  and 
Mechanical  Engineering 
Tucson,  Arizona  85721 

Dr.  S.  J.  Veuves 
Carnegie-Mellon  University 
Department  of  Civil  Engineering 
Schenley  Park 

Pittsburgh,  Pennsylvania  15213 

Dr.  Bonald  L.  Huston 
Department  of  Engineering  Analysis 
University  of  Cincinnati 
Cincinnati,  Ohio  45221 

Professor  6.  C.  M.  Sih 
Lehigh  University 
Institute  of  Fracture  and 
Solid  Mechanics 
Bethlehem,  Pennsylvania  18015 

Professor  Albert  S.  Kobayashi 
University  of  Washington 
Department  of  Mechanical  Engineering 
Seattle,  Washington  98105 

Professor  Daniel  Frederick 
Virginia  Polytechnic  Institute  and 
State  University 

Department  of  Engineering  Mechanics 
Blacksburg,  Virginia  24061 
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Universities  (Con't) 

Professor  A.  C.  Eringen 
Princeton  University 
Department  of  Aerospace  sad 
Mechanical  Sciences 
Princeton,  Hew  Jersey  08540 

Professor  E.  H.  Lee 

Stanford  University 

Division  of  Engineering  Mechanics 

Stanford,  California  94305 

Professor  Albert  I.  King 
Wayne  State  University 
Biomechanics  Kesearch  Center 
Detroit,  Michigan  48202 

Dr.  V.  L  Hodgson 
Wayne  state  University 
School  of  Medicine 
Detroit,  Michigan  48202 

Dean  B.  A.  Bo ley 
Northwestern  University 
Department  of  .Civil  Engineering 
Evanston,  Illinois  60201 


Professor  H.  M.  Nevmark 
University  of  Illinois 
Department  of  Civil  Engineering 
Urbane,  Illinois  61803 


Professor  E.  Beissner 
University  of  California,  San  Diego 
Department  of  Applied  Mechanics 
La  Jolla,  California  92037 

Professor  William  A.  Hash 
University  of  Massachusetts 
Department  of  Mechanics  and 
Aerospace  Engineering 
Amherst,  Massachusetts  01002 


Professor  P.  C.  Hodge,  Jr. 
University  of  Minnesota 
Department  of  Aerospace  Engineering 
and  Mechanics 

Minneapolis,  Minnesota  55455 

Dr.  D.  C.  Drucker 
University  of  Illinois 
Dean  of  Engineering 
Urbane,  Illinois  61801 
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Universit  1m  (Con't) 

Professor  G.  Herraenn 
Stanford  University 
Department  of  Applied  Machanies 
Stanford,  California  94305 

Profasaor  J.  D.  Achanbacb 
Northwest  Uni varsity 
Department  of  Civil  Engineering 
Evanston,  Illinois  60201 

Profasaor  S.  B.  Dong 
Univarsity  of  California 
Department  of  Machanies 
Los  Angolas,  California  90024 

Profasaor  Burt  Paul 
Univarsity  of  Pannsylvania 
Towns  School  of  Civil  and 
Mechanical  Enginaaring 
Philadelphia,  Pennsylvania  19104 

Professor  H.  W.  Liu 
Syracuse  Univarsity 
Department  of  Chemical  Engineering 
and  Metallurgy 
Syracuse,  Maw  York  13210 

Professor  S.  Bodnar 
Technion  BAD  Foundation 
Haifa,  Israel 

Professor  Werner  Goldsmith 
Univarsity  of  California 
Dapartaant  of  Mechanical  Engineering 
Berkeley,  California  94720 

Professor  R.  S.  Rivlin 
Lehigh  Univarsity 
Canter  for  the  Application 
of  Mathematics 

Bethlehem,  Pannsylvania  18013 

Profasaor  P.  A.  Coxsarelli 
State  Univarsity  of  New  Tork  at 
Buffalo 

Division  of  Interdisciplinary  Studies 
Ran  Parker  Engineering  Building 
Chemistry  Road 
Buffalo,  Mew  Tork  14214 


Universities  (Con't) 

Professor  Joseph  L.  Rose 
Dr axel  University 

Department  of  Mechanical  Engineering 
and  Mechanics 

Philadelphia,  Pennsylvania  19104 

Professor  8.  R.  Donaldson 
University  of  Maryland 
Aerospace  Engineering  Department 
College  Park,  Maryland  20742 

Professor  Joseph  A.  Clark 
Catholic  University  of  America 
Department  of  Mechanical  Engineering 
Washington,  D.C.  20064 

Dr.  Samuel  B.  Batdorf 
University  of  California 
School  of  Engineering 
and  Applied  Science 
Los  Angeles,  California  90024 

Professor  Isaac  Fried 
Boston  University 
Depertment  of  Mathematics 
Boston,  Massachusetts  02215 

Professor  E.  Rrempl 
Rensselaer  Polytechnic  Institute 
Division  of  Engineering 
Engineering  Mechanics 
Troy,  Hew  York  12181 

Dr.  Jack  R.  Pinson 
University  of  Delawere 
Department  of  Mechanical  and  Aerospace 
Engineering  and  che  Center  for 
Composite  Materials 
Newark,  Delaware  19711 

Dr.  J.  Duffy 
Brown  University 
Division  of  Engineering 
Providence,  Rhode  Island  02912 

Dr.  J.  L.  Swedlow 
Camegie-Mellon  University 
Department  of  Mechanical  Engineering 
Pittsburgh,  Pennsylvania  15213 
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Universit ias  (Con't) 

Or.  V.  K.  Varadan 

Ohio  Scata  University  Raaearch  Foundac icm 
Department  of  Engineering  Mechanics 
Columbus,  Ohio  43210 

Or.  Z.  Hash in  .  ■ 

University  of  Pennsylvania 
OapartaanC  of  Metallurgy  and 
Materials  Science 
College  of  Engineering  and 
Applied  Science 

Philadelphia,  Pennsylvania  19104 

Dr.  Jackson  C.  S.  Tang 
University  of  Maryland 
Oepartaent  of  Mechanical  Engineering 
College  Park,  Maryland  20742 

Professor  T.  T.  Chang 
University  of  Akron 
Oepartaent  of  Civil  Engineering 
Akron,  Ohio  44325 

Professor  Charles  W.  Bert 
University  of  Oklahoma 
School  of  Aerospace,  Mechanical, 
and  Nuclear  Engineering 
Noraan,  Oklahoma  73019 

Professor  Satya  N.  Atluri 
Georgia  Institute  of  Technology 
School  of  Engineering  and 
Mechanics 

Atlanta,  Georgia  30332 

Professor  Graham  P.  Carey 
University  of  Texas  at  Austin 
Oepartaent  of  Aerospace  Engineering 
and  Engineering  Mechanics 
Austin,  Texas  78712 

Dr.  S.  S.  Wang 
University  of  Illinois 
Oepartaent  of  Theoretical  and 
Applied  Mechanics 


Universities  (Con't) 

Professor  V.  H.  Neubert 
Pennsylvania  State  University 
Department  of  Engineering  Science 
and  Mechanics 

University  Park,  Pennsylvania  16802 

Professor  A.  W.  Leissa 
Ohio  State  University 
Department  of  Engineering  Mechanics 
Columbus,  Ohio  43212 

Professor  C.  A.  Brebbia 
University  of  California,  Irvine 
Department  of  Civil  Engineering 
School  of  Engineering 
Irvine,  California  92717 

Dr.  Ceorge  T.  Hahn 
Vanderbilt  University 
Mechanical  Engineering  and 
Materials  Science 
Nashville,  Tennessee  37235 

Dean  Richard  H.  Gallagher 
University  of  Arizona 
College  of  Engineering 
Tucson,  Arizona  85721 

Professor  E.  7.  Rybicki 
The  University  of  T'tlsa 
Department  of  Mechanical  Engineering 
Tulsa,  Oklahoma  74104 

Dr.  R.  Haftka 

Illinois  Institute  of  Technology 
Department  of  Mechanics  and  Mechanical 
and  Aerospace  Engineering 
Chicago,  Illinois  60616 

Professor  J.  G.  de  Oliveira 
Massachusetts  Institute  of  Technology 
Department  of  Ocean  Engineering 
77  Massachusetts  Avenue 
Cambridge,  Massachusetts  02139 


Urbane,  Illinois  61801  Dr.  Bernard  tf.  Shaffer 

Polytechnic  Institute  of  New  Tork 
Professor  J.  7.  Abel  Route  110 

Cornell  University  Paraingdale,  New  Tork  11735 

Oepartaent  of  Theoretical 
and  Applied  Mechanics 
Ithaca,  New  Tork  14853 
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Industry  and  Research  Institute 


Or.  Norman  Hobbs 
Kaman  AviDyue 
Division  of  Kama n 
Sciancea  Corporation 
Burlington,  Maaaaehuaatts  01803*  - 

Argonna  National  Laboratory 
Library  Sarvices  Department 
9700  South  Caaa  Avenue 
Argonna,  Illinois  60440 

Dr.  N.  C.  Juager 
Cambridge  Acoustical  Aasociates 
34  Rindge  Avenue  Extension 
Cambridge,  Massachusetts  02140 

Mr.  J.  H.  Torrance 
Ceneral  Dynamics  Corporation 
Electric  Boat  Division 
Croton,  Connecticut  06340 

Dr.  J.  E.  Greenspon 
J.  G.  Engineering  Research  Associates 
3831  Menlo  Drive 
Baltimore,  Maryland  21215 

Newport  News  Shipbuilding  and 
Dry  Dock  Company 
Library 

Newport  News,  Virginia  23607 

Dr.  W.  F.  Boxich 

McDonnell  Douglas  Corporation 

5301  Boise  Avenue 

Huntington  Beach,  California  92647 

Dr.  H.  N.  Abramson 
Southwest  Research  Institute 
8300  Culebra  Road 
Sea  Antonio,  Texas  78284 

Dr.  R.  C.  DeHart 
Southwest  Research  Institute 
8500  Culebra  Road 
San  Antonio,  Texas  78284 

Dr.  M.  L.  Baron 
Weidlinger  Associates 
110  East  39th  Straec 
New  Tork,  New  York  10022 


Industry  and  Research  Institutes  (Con't) 


Dr.  T.  L.  Geers 

Lockheed  Missiles  and  Space  Company 

3251  Hanover  Street 

Palo  Alto,  California  94304 

Mr.  William  Cayvood 
Applied  Physics  Laboratory 
Johns  Hopkins  Road 
Laurel,  Maryland  20810 

Dr.  Robert  E.  Dunham 
Paeifica  Technology 
P.0.  Box  148 

Del  Mar,  California  92014 

Dr.  M.  P.  Kaaninen 
Battalia  Columbus  Laboratories 
SOS  King  Avenue 
Columbus,  Ohio  43201 

Dr.  A.  A.  Hochrein 
Daedalean  Associates,  Inc. 
Springlake  Research  Road 
15110  Frederick  Road 
Woodbine,  Maryland  21797 

Dr.  James  W.  Jones 
Swanson  Service  Corporation 
P.0.  Box  5415 

Huntington  Beach,  California  92646 

Dr.  Robert  E.  Nickell 
Applied  Science  and  Technology 
3344  North  Torrey  Pines  Court 
Suite  220 

La  Jolla,  California  92037 

Dr.  Kevin  Thomas 
Wcstinghouse  Electric  Corp. 

Advanced  Reactors  Division 
P.  0.  Box  158 

Madison,  Pennsylvania  15663 

Dr.  H.  D.  Hibbitt 
Hibbitt  &  Karla son,  Inc. 

132  George  M.  Cohan  Boulevard 
Providence,  Rhode  Island  02903 

Dr.  R.  D.  Mindlin 
89  Deer  Hill  Drive 
Ridgefield,  Connecticut  06877 
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Industry  and  Research  Institute » 

Or.  Richard  E.  Dana 
Mega  Engineering 
11961  Tech  Road 

Silver  Spring,  Maryland  20904 

Mr.  C.  M.  Stanley 
Lockheed  Palo  Alto  Research 
Laboratory 
3251  Hanover  Street 
Palo  Alto,  California  94304 

Mr.  R.  L.  Cloud 

Robert  L.  Cloud  Associates,  Inc. 
2972  Adeline  Street 
Berkeley,  California  94703 


(Con’t) 


DATE 

ILME 


